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FOREV/ORD 


This  second  Technical  r>ocumcntary  Report  presents  a  technical  disc  ission 
and  review  of  program  activity  and  progress  during  the  period  16  September 
through  15  December  1964. 
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(This  abstract  is  unclassified) 


ABSTRACT 


The  results  of  an  investigation  of  viscoelastic  and  failure  properties  of  highly 
filled  FbAA  and  FBAN  propel  Irtiits  as  a  function  ol  ooiiuo  leading  are  lepoiiei 
Failure  surface  study  results  are  reported,  and  the  results  of  a  limited  study 
of  the  relationship  between  crack  propagation  velocity  and  propellant  physical 
characteristics  are  discussed.  Propellant  dynamic  shear  and  bulk  properties 
were  investigated  with  small  deformation  piezoelectric  drvicer.  An  experi¬ 
mental  investigation  of  propellant  thermomechanical  response  to  sustained 
cyclic  inertial  loading  was  completed,  and  the  results  in  agreement  with  theoi 
are  presented.  Also  di.scussed  are  experbriental  investigations  of  transient 
thermoviscoelastic  response  of  propellant  under  c  "instant  cyclic  strain  ampli¬ 
tude  and  inertial  loading. 
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Section  1 

INTRODUCTION  ANT)  SUMMARY 


The  Structural  Integrity  Department  of  Lockheed  Propulsion  Company  is 
conducting  a  theoretical  and  experimental  research  program  to  correlate 
the  mcclianical  behavior  and  failure  properties  of  solid  propellant  with  the 
basic  characteristics  of  the  propellant  and  its  constituents.  This  research 
is  supported  by  the  Air  Force  Rocket  Propulsion  Laboratory,  Research  and 
Technology  Division,  Air  Force  Systems  Command,  Edwards,  California, 
under  Contract  Ko.  AF  04(6 1 1) -9953.  The  program  effort  is  concentrated 
on  the  study  of  propellant  dynamic  physical  properties,  structural  failure 
phenomena,  and  chemical-physical  effects. 

In  this  report,  program  activities  and  results  for  the  period  15  September 
through  15  December  1964  are  described  and  discussed.  A  section-by¬ 
section  summary  is  presented  below. 

Section  2  —  Propellant  Formulation  StufL/,  lligh  .Solids  Loading.  Fc'iilure 
proj  erty  investigations  in  mulTTa'xTal  stre.s.s  states  were  completed  for  the 
.second  group  of  PBAA  and  PLAN  propell.ints  which  include  oxidizer  modality 
and  burning  rate  catalyst  variations.  Little  effect  on  failure  properties  was 
observed  for  the  change  from  optimized  trioxidizer  to  optimized  bioxidizer 
particle  size  distribution  as  anticipated  from  theoretical  considerations.  Some 
sensitivity  of  the  failure  properties  of  the  PBAA  propellants  to  the  liquid  burn 
rate  catalyst  was  observed.  Stress  relaxation  characterization  and  small  strain 
nonlinear  behavior  of  the  first  group  of  formulation  variations  are  discussed 
along  with  prelimin  ry  correlation  of  the  relaxation  data  with  measured  small 
deformation  dynamic  response. 

Section  3  —  Fracture  M «> c iia n i c s .  Work  during  the  quarter  was  concentrated 
on  lalivrre  .surface  cTefTrutTon  in  orthogonal  stress  space  and  studies  of  crack 
propagation  in  biaxial  strip  test  specimens.  Studies  of  test  specimens  for 
obtaining  failure  data  in  three-dimensional  stress  states  have  continued.  Re¬ 
sults  of  the  failure  surface  studies  are  presented.  A  preliminary  estimate  has 
been  made  of  the  shape  of  the  failure  surface  in  stress  space  for  an  84  weight- 
percent  solids-loadcd  PBAN  propellant.  A  maximum  tensile  stress  criterion 
appears  consisti.nt  with  the  limited  data  for  t)  e  tensile  octant;  an  internal 
friction  criterion  appears  generically  applicable  to  the  octants  where  tensions 
and  compressions  arc  combined.  Superposed  hydrostatic  pressure  on  uniaxial 
tension  caused  only  small  changes  in  strength  and  elongation. 

The  results  of  the  limited  .studies  and  experiments  that  have  directed  attention 
to  relationships  between  crack  propagation  velocity  and  propellant  physical 
characteristics  are  described.  Good  agreement  was  found  between  measured 
and  predicted  crack  exten.sions. 
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Section  4  —  Dynamic  Properties.  Experimental  investigations  of  propellant 
<i y n a m i c *s h e arlinirFulK  properties  measured  Vk'ith  the  small  deformation 
piezoelectric  devices  previously  described  (Ret.  1  )  are  discussed.  Shear 
moduli  of  PBAA  and  PBAN  propellant  and  cured  binder  are  presented  and 
compared  with  linear  thermal  expansion  measurements  of  glass  transition 
temperature  obtained  for  similar  propellants.  Analysis  techniques  and  cali¬ 
bration  results  are  discussed  for  the  dynamic  bulk  modulus  device. 

An  experil  umtal  investigation  of  propellant  thermomechanical  response  to 
sustained  i.yclic  inertial  loading  was  completed.  Experimental  results  are 
presented  and  it  is  noted  that  they  are  in  qualitative  agreement  with  theoreti¬ 
cal  predictions.  In  particular,  jump  instability  was  experimentally  observed 
to  be  due  to  regenerative  thermal  and  mechanical  coupling.  Also,  an  experi¬ 
mental  investigation  of  transient  thermoviscoelastic  response  for  constant 
cyclic  strain  amplitude  and  inertial  loading  was  completed.  Experimental 
results  presented  include  cyclic  degradation  observations. 
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Section  2 


PROPELLANT  FORMULATION  STUDY 


HIGH  SOLIDS  LOADIl 


2.1  INTRODUCTION 


Increasing  use  is  being  made  of  very  highly  loaded  solid  propellants  in  prac¬ 
tical  motor  applications.  As  the  solids  loading  is  increased,  the  physical 
properties  change  and  the  propellants  become  more  nonlinear  in  their  visco¬ 
elastic  behavior.  For  these  reasons,  research  at  Lockheed  Propulsion 
Company  in  the  areas  of  dynamic  response  and  failure  laechanisms  is 
presently  being  concentrated  on  ballistically  feasible  propellant  formulations 
approaching  the  maximum  theoretical  limit  of  solids  loading.  Primary 
effort  during  the  second  quarterly  reporting  period  has  been  directed  towards 
completion  of  an  experimental  parametric  study  of  a  series  of  highly  loaded 
PBAA  and  PBAN  propellants. 

Although  analysis  of  the  experimental  results  are  so  far  incomplete,  several 
conclusions  concerning  the  observed  effects  of  high  rsolids  loading  on  propel¬ 
lant  physical  properties  can  be  made  at  this  time,  A  systematic  degradation 
of  the  propellant  failure  strain  capability  and  an  accompanying  increase  in 
failure  stress  are  observed  as  the  solids  loading  is  increased,  resulting  in 
a  net  increase  in  the  mechanical  energy  required  to  break  the  specimens.  In 
addition,  the  viscoelastic  properties  at  small  strain.?  are  observed  to  mirror 
the  failure  behavior  with  a  general  tendency  toward  higher  moduli  and  increas¬ 
ing  nonlinear  behavior  as  the  solids  loading  increases. 

Experimental  results  for  the  burning  rate  catalyzed  and  bimodal  oxidizer 
propellant  formulations  are  presented  in  the  following  sections.  The  results 
show  that  changing  the  oxidizer  modality  from  a  trioxidizer  to  a  hioxidizer 
system  has  little  effect  on  the  failure  properties  when  the  particle  size  dis¬ 
tributions  for  the  tw'o  modal  systems  are  optimized  as  discussed  in  the 
previous  quarterly  report. 

In  general,  burning  rate  catalysts  were  found  to  ha\  e  little  effect  on  failure 
properties  with  the  single  exception  that  PBAA  propellants  were  found  to 
be  softened  substantially  by  the  addition  of  the  liquid  burning  rate  catalyst. 
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2.2  TEST  PLAN 


Propellants  have  been  prepared  consisting  of  combinations  of  four  propellant 
formulation  variables.  These  arc  (1)  total  weight  percent  solids  loading, 

(2)  binder  type,  (3)  oxidizer  particle  size  distribution,  and  (4)  burning  rate 
catalysts.  The  total  solids  loading  was  varied  between  84  and  90  percent. 
Polybutadiene  acrylic  acid  (PBAA)  copolymer  and  polybutadiene  acrylonitrile 
acrylic  acid  (PBAN)  terpolymer  binder.^  v.’ere  investigated.  Two  particle 
size  distribution  systems  were  used--a  bimodal  oxidizer  system  and  a  tri- 
modal  oxidizer  system,  both  optimized  for  minimum  processing  viscosity. 

In  addition,  an  assessment  of  the  relative  effects  of  burning  rate  catalysts, 
both  liquid  and  solid,  on  the  physical  properties  and  burning  rate  was  planned 
on  a  spot-check  basis.  Formulations  which  have  been  cast  and  prepared  for 
testing  are  listed  in  Table  2-  1.  Formulations  will  be  referred  to  by  the  code 
symbols  listed  in  Table  2-1  throughout  this  report.  All  the  propellants  ^re 
aluminized.  The  solids  particle  size  distribution  is  identical  for  each  tri¬ 
oxidizer  System.  Similarly,  each  bioxidizer  system  contains  an  identical 
solids  system,  with  the  exception  of  the  two  systems  catalyzed  with  solid 
catalysts,  84NBCS  and  88NBCS.  For  these  systems,  a  small  portion  of  the 
fine  oxidizer  is  replaced  with  an  equivalent  amount  of  burning  rate  catalyst, 
PBAA  binder  is  cured  at  a  l.l/l  imine-to-carboxyl  ratio  level  with  a 
trifunctional- to- difunctional  curative  ratio  of  40/60.  PBAN  binder  is  cured 
at  a  1.3/1  epoxy-to-carboxyl  ratio.  Both  PBAA  and  PBAN  propellant  binders 
are  plasticized.  Specific  information  regarding  the  formulation  of  each  pro¬ 
pellant  was  presented  in  Appendix  II  (Confidential)  of  LPC  Report  667-Q-l 
(AFRPL-TR-64-  148). 

Two  5500-gram  batches  of  each  formulation  were  mixed  simultaneously  in 
similar  laboratory  mixing  facilities  to  provide  sufficient  propellant  for  the 
various  specimens  and  to  eliminate  batch-to-batch  variations  as  much  as 
possible.  The  propellants  were  mixed  and  cast  under  vacuum. 

The  test  specimens  used  in  physical  property  testing  were  the  in  situ  bonded 
tab-end  type,  sawed  and  milled  from  bulk  casting  of  propellant.  Prior  to 
testing, all  specimens  were  dimensionally  inspected  and  stored  in  a  70®F 
desiccated  environment. 

Experimental  data  for  four  stress  states  are  being  obtained  for  the  propellant 
formulaticns  in  Table  2-1  .  Uniaxial  tension,  strip  biaxial,  diametral  com¬ 
pression,  and  double-lap  shear  tests  are  being  performed,  using  a  constant 
rate  of  deformation  for  failure  characterization. 

The  essential  characteristics  of  the  specimens  and  the  definitions  of  failure 
are  discussed  in  the  following  subsections. 
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2.2.1  Uniaxia  1  Tension’ 


Tab-end  specimens  consisting  of  rectangular  bars  (4  by  O.J  by  C.B-incbcs) 
with  wooden  tab-endn  (bonded  in  situ  to  propellant  during  cure)  are  being  used 
exclusively  for  the  tensile  tests.  The  stre-ss-strain  condition  identified  as 
failure  i.®.  defined  by  <  the  so-calkd  strain  at  maximuni  stress,  as  pre¬ 
sented  in  the  ICRPG  Solid  I^ropellant  h’echanita  1  Behavior  Manua'  (Ilcf.  2]. 
This  strain  fevcT  i.s  usually  found  to  bT~aTFout~Z7*3''or  tlic7  .stFaTir"icVtT  where  the 
load-time  record  passes  through  a  maximum,  and  may  be  associated  with 
the  development  of  extensive  dewetting. 


2.2,2  Strip  Biaxial' 


Strip  biaxial  .specimens  consist  of  thin  plates  of  propellant  with  a  width  much 
greater  than  their  hcigVit.  The  specimens  are  prepared  by  casting  piopellant 
into  lined  wooden  boxes  of  appropriate  size  and  subsequently  machining 
specimens  from  the  box,  retaining  the  w'ooden  tabs  along  the  long  sides  of 
the  specimens.  They  are  tested  at  a  constant  elongation  rate  applied  in  a 
direction  normal  to  the  clamped  edges.  In  most  instances,  strip  specimens 
with  a  large  width-to-height  ratio  fail  by  cracking  in  the  central  region  of 
the  strip  where  the  stress  field  iw  reasonably  uniform.  As  in  the  uniaxial 
tensile  tests,  a  maximum  in  the  load-time  curve  sometimes  occurs  prior  to 
visual  crack  initiation.  Thus,  for  the  sake  of  uniformity,  the  stress- strain 
condition  identified  as  failure  in  this  test  is  defined  in  the  same  manner  as 
for  the  unia.xial  tensile  tests.  The  size  of  the  biaxial  strip  specimens  used 
is  7,5  by  15  by  0,125  inches,  and  the  sheets  arc  filleted  from  the  0. 12 5 -inch 
thickness  to  a  3/8-inch  thick  wooden  tab  strip. 

2.2,3  Diametral  Compression’ 


The  diametral  compression  test  specimens  arc  1.5-inch  diameter  discs, 
l/?.-inch  thick,  prepared  by  machining  cast  cylinde.-s  of  pronellp»nt.  The 
failure  of  these  specimens  is  first  evidenced  by  the  appearance  of  fissures 
at  the  center,  on  the  surface  of  the  specimens,  usually  aligned  at  an  angle 
to  the  principal  stress  axes  so  as  to  suggest  shear  failure.  Although  it  is 
possible  that  the  initiation  of  failure  takes  place  inside  the  sample  somewhat 
earlier  than  is  shown  by  the  surface  crack,  the  start  of  the  visible  fissuring 
at  the  specimen  surface  defines  failure. 


’These  test  specimens  are  described  in  detail  in  Reference  2J, 
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2.2.4  Double 5>bear 


Parallelogram  plan -form  shear 
lant  failure  properties  undei 
supports  constrained  to  pre^ 
the  center  svippcrt  rc.'ralis 
sequent  failure  of  the  prop 
Shear  failure  surfaces  are 
perpendicular  tc  the  free  surfaces  of  the  propellant.  Visual  observation  of 
fracture  is  difficult,  br.t  it  is  usually  observed  nearly  at  the  onserved  peak 
output  force  in  a  constant  shear  strain  rate  test.  For  simplicity,  therefore, 
failure  is  defined  in  this  test  as  the  peak  in  the  force-time  record. 
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2.3  FAILURE  PROPERTIES 


Multiaxial  failure  characteristics  of  the  10  fortriulations  shown  unshaded  in 
Table  2-  '.  have  been  investigated.  The  results  show  that  the  absolute 
magnitudes  of  the  failure  stress  and  strain  again  vary  with  stre5ja  axiality 
as  anticipated.  There  is  also  a  consistent  change  in  failure  prop>crties  with 
formulation  variation  for  all  of  the  test  modc.s.  Thus,  hccau;'c  of  the  largo 
amount  of  d-ata,  unnecessary  repetitio.n  will  ViC  jivoidod  by  prr  r. r-r.ting  in 
graphical  form  only  the  strip  biaxial  test  mode  rcsu''t3.  railurc  data  fen* 
the  remaining  three  tc.st  modes  are  presented  in  tabular  form  in  Tables  2-2  , 
2-3  ,  and  2-4  . 

Figure  2-  I  presents  comparative  biaxial  strip  failure  properties  for  the  84 
weight  •  per  rent  loaded,  bioxidizer,  PBAN  catalysed  and  uncalalyzed  formu¬ 
lations  84NBC3  and  84  NB,  respectiv-ely.  The  solid  burn  rate  cnt.alyzing 
agent  (less  than  two  percent  by  weight  of  the  total  formulation)  replaced  an 
equivalent  weight  of  fine -ground  ammonium  perchlorate  to  maintain  a  con¬ 
stant  solids  loading.  Inspection  of  the  data  showed  a  small  amount  of  pro¬ 
pellant  stiffening  due  to  the  solid  catalyst  resulting  in  slightly  higher  failure 
stresses  and  lower  failure  strains.  This  effect  is  confirmed  by  the  diametral 
compression  and  lao-shear  test  data. 

Similar  PBAN,  bioxidizer  formulations  at  88  percent  solids  loading,  with 
various  burning  rate  catalysts,  are  compared  in  Figure  2-2.  No  significant 
change  in  failure  'Toperfies  is  apparent  from  these  data  as  a  result  of  using 
the  burning  rate  catalysts.  The  liquid  catalyst  amounts  to  less  than  two 
percent  by  weight  of  the  total  formulation  and  replaces  an  cqu.al  quantity  of 
plastici'zer  so  as  to  keep  the  total  plasticizer -to -polymer  ratio  constant. 
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Table  2-2 

UNIAXIAL  FAILURE  STRESS  AND  STRAIN 
AT  ICRPG  (1963)  FAILURE  POINT,  Rate  5%/min 


Failure  Stresn  (psi) 
TSHurb  Strain  (percent) 


Propellant 

Formulation 

Temperature  (®F) 

140 

83 

10 

-50 

-85 

84  RB 

55 

79 

212 

38  3 

635 

5.3 

6.0 

7,6 

5.1 

4.81 

84  NB 

53 

76 

138 

421 

850 

19 

22.1 

12.5 

6.2 

2.95 

88  RB 

79 

124 

294 

-- 

-- 

4.1 

4.3 

5.2 

-- 

-- 

88  NB 

72 

96 

190 

557 

773 

13 

13.1 

12.1 

5.3 

1.71 

84  NBCS 

69 

91 

185 

485 

802 

18 

18.6 

12.6 

7.1 

2.38 

88  NBCS 

88 

108 

196 

538 

882 

14 

13.8 

11.7 

5.6 

2.25 

88  RC 

31 

54 

167 

326 

457 

8.2 

8.4 

8.2 

3.08 

2.37 

88  NC 

66 

87 

181 

530 

662 

WBm 

9.2 

9.2 

3.20 

2.00 

88  RBC 

40 

68 

164 

326 

446 

14.2 

12.7 

9.7 

4.72 

2.99 

1 

88  NBC 

81 

107 

198 

543 

-- 

10.7 

9.7 

6.7 

3.43 

-- 
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Table  2-3 

DIAMETRAL  COMPRESSION  HORIZONTAL  FAILURE  STRESS  AM 
AT  CENTER,  Rate  6.3%/min  VERTICAL  DIAMETER  CHAN 


Failure  Stress  (psi) 
Failure  Strain  (percent) 


Propellant 

Temperature  {°F) 

Formulation 

140 

74 

10 

84  RB 

52.3 

102 

2  35 

20.3 

27.2 

31.2 

84  NB 

43.6 

64.8 

211 

37.1 

53.6 

56.9 

88  RB 

63.9 

127.6 

266 

13.1 

18.0 

20.0 

68  NB 

54.0 

81.7 

163 

27.0 

30.1 

36.2 

84  NBCS 

53.2 

91.7 

213 

35.7 

44.3 

,53.7 

88  NBCS 

69.5 

100 

170 

24.2 

28.2 

26.5 

88  RC 

30.6 

102 

196 

28.5 

30.1 

29.1 

88  NC 

52.3 

77.9 

170 

25.7 

22.0 

28.0 

88  Rnc 

38.7 

89.6 

169 

33.9 

31.3 

32.9 

88  NRC 

71.4 

112 

240 

1 _ 

19.5 

25.2 

26.6 
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Table  2-4 

DOUBl  E-LAP  SHEAR  FAILURE  STRESS  AND  STRAIN. 

Rate  6 .7%/min 


Failure  Stress  (psi) 
Failure  Strain  (percent) 


Propcllc-nt 

Formulation 

Temperature  (°F) 

140 

82 

84  RB 

50,2 

80.9 

23.9 

25.7 

84  NB 

40.5 

62.1 

41.5 

4  (.8 

88  RB 

54.0 

89.2 

15.1 

18.0 

88  NB 

60.9 

f  0.1 

31.8 

34,8 

84  NBCS 

38.2 

66.5 

33,9 

37.4 

88  NBCS 

64.9 

104 

31.1 

34,1 

88  RC 

29,2 

75.1 

34.1 

34.1 

88  NC 

49.6 

88.6 

21.3 

31.2 

88  RBC 

25.9 

68.5 

39.4 

37.5 

88  NBC 

65.8 

113 

25.1 

32.2 
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Fipuic  ?- I  Biaxial  Strij'  Failure  Stress  anrl  Strain  at  ICRPG 
(1963)  Failure  Point  versus  Tojnperature,  PBAN, 
34^  L  o  a '  1 :  n ,  13  i  o i  i  j  X,  e  r 


2-9 


LfteifUBvn 


AVERAGE  STRAIN  «Pei»CtNT)  AVERAGE  STRESS  (PSi» 


1*  i{,urc  2-2  Biaxial  Strip  Failure  Stress  and  Strain  at  ICRPG 
(1963)  Failure  Point  versus  Tejr.pc raturo,  PBAN, 
88)S  Lead  ing,  Bioxidizer 
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In  Figure  2-3,  the  effects  of  solids  loading  on  the  PBAA  and  PBAN  bioxicliuer 
formulations  are  compared  at  the  84  and  88  weight -percent  solids  loading 
levels  .  The  failure  characteristics  are  similar  to  the  trioxidiacr  formula¬ 
tions  in  that  increasing  solids  loading  results  in  increased  failure  stresfs 
and  reduced  failure  strains. 

Brittle  fracture  was  obtained  with  the  PBAN  propellants  when  tested  at 
-50  and  -85®F.  Again  it  was  found  that  the  PBAA  propellants  have  unusually 
low  failure  strains  and  high  failure  stresses  as  noted  for  the  uncatalyxed 
trioxidizer  PBAA  formulations. 

A  comparison  between  the  uncatalyzed  trioxidizer  and  bioxidizer  formula¬ 
tions  is  presented  in  Figures  2-4  and  2-5  for  the  PBAN  and  PBAA  propel¬ 
lants,  respectively.  In  general,  the  bioxidizer  formulations  show  slightly 
superior  strain  and  stress  capability  for  both  binder  types  and  solids  load¬ 
ing  le-.-elc  .  However,  this  situation  is  reversed  at  higher  temperatures 
with  the  PBAA  formulations. 

Figures  2-6  and  2-7  illustrate  the  effect  of  the  liquid  burning  rate  catalyst 
on  the  failure  properties  of  the  bioxidizer  and  trioxid-zer  formulations  at 
88  percent  loading  for  F'BAA  and  PBAN  propellants,  respectively.  The 
dat?  for  the  PBAA  binder  show  a  significant  softening  of  the  propellant  when 
the  liquid  burning  rate  catalyst  is  added.  Failure  strain  capability  is  in¬ 
creased  by  approximately  a  factor  of  two  at  the  higher  temperatures  and 
the  stresses  are  correspondingly  lower.  Inspection  of  Figure  2-7  reveals 
no  such  effect  for  the  PBAN  binder  propellants  and  in  fact  a  slight  stiffening 
is  observed  due  to  the  liquid  catalyst.  The  data  therefore  show  a  preferent¬ 
ial  interaction  between  the  liquid  burning  rate  catalyst  and  the  PBAA  binder, 
tending  to  reduce  the  effective  network  crosslinking  or  otherwise  loosen  the 
binder  structure.  The  fact  that  the  biggest  effect  is  obtained  at  higher 
temperatures  suggests  that  binder  cros slinki.ag  structure  (or  effective 
chain  length)  is  modii.ed  by  the  catalyst,  resulting  in  a  lower  modulus  in 
the  temperature  region  where  long  range  interactions  are  most  important. 


2-11 


AVERAGE  STRAIN  l««C«MTI  AVERAGE  STRESS  1^1 


AFRPL-TR-65-20 


667.Q-2 


Fif  ure  2-3  Biaxiai  Sirip  Failure  Stress  and  Strain  at  ICRPG 
(1963)  Failure  Point  vers  IS  Temperature, 
Bioxidi.jer,  Ui.'.“a»^aly7.ed 
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Fif'iii-f  ^'-5  liiaxitil  Slrip  f'ailure  Stress  and  Strain  at  ICRPG 
(1963)  Failare  Point  versa  i  Temperature,  PBAA, 
U  nr  a  IS  ly/.i'f' 
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Figure  2-6  Biaxial  Strip  Failure  Stress  a, id  Strain  at  ICRPG 
(1963)  Failure  Point  versus  Temperature,  PBAA, 
88^  Loading 
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Figure  2-V  Biaxi.'J  Strsp  )"-.rlnre  Sliess  and  otrain  at  ICRPG 
(1963)  Failure  Point  versus  Temperature,  PBAM, 
Loading 
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for  tlie  remainirig  high  soIiJs  loaded  propcllnnls  eLowed  that  the  relative 
magnitvide  of  tliia  particular  effect  increase.^  as  the  eolids  loading  is  increasec! 
and  as  the  binder  modulus  decreases  at  a  constant  solids  londing. 

The  study  of  the  tensile  relaxation  response  of  the  first  group  of  highly  loaded 
propellant  formulations  has  been  continued  and  extended  to  temperatures 
ranging  from  -24°  F  to  +174°F.  In  general,  relaxation  response  was  found 
to  be  a  nonlinear  function  of  strain  level  in  the  1  to  5  percent  strain  range. 

In  addition,  time -temperature  superposition  of  the  data  was  found  to  be  poor. 
However,  mathematical  inversion  of  the  reduced  relaxation  data  to  complex 
tensile  modvilus  form  was  performed  and  the  results  were  compared  with 
frequency-temperature  superposed,  experimentally  measured  dynamic  shear 
data  for  the  same  propellant.s.  The  dynamic  shear  tests  were  performed  at 
very  low  strains  of  the  order  of  6  x  10“*  percent.  In  general,  results  which 
have  been  analyzed  thus  far  show  good  agreement  betw'ecn  the  loss  moduli 
components  from  the  tensile  relaxation  data  and  the  measured  dynamic  shear 
data.  However,  the  storage  component  curves  do  not  show  good  agreement. 
The  storage  moduli  calculated  from  the  tensile  i'claxation  tests  are  generally 
larger  than  the  values  calculated  from  the  sni.iU  deformetien  shear  data.  I'lie 
discrepancy  increases  at  lower  frequencies  or  at  higher  temperatures,  result¬ 
ing  in  a  slope  micmatch  between  the  curves  of  log  modulus  vc.r.sus  log  reduced 
frequency  in  the  overlap  region. 

Data  reduction  and  analysis  are  still  in  progress  and  detailed  results  will  be 
given  in  subsequent  reports. 
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FRACTURE  MECflANZCS 


3.1  GENERAL 


Work  during  the  quarter  was  concentrated  on  failure  surface  definition  in 
orthogonal  stress  space  and  studies  of  crack  propagation  in  biaxial  strip  test 
specimens.  Studies  of  test  specimens  for  obtaining  failure  data  in  three 
dimensional  stress  states  have  continued.  Section  3,2  presents  the  results 
of  the  failure  surface  studies.  The  inform?tion  in  this  section  is  essentially 
identical  to  that  contained  in  AIAA  preprint  paper  No.  65-157.  A  preliminary 
estimate  has  been  made  of  tlie  shape  of  the  failure  surface  in  stress  space 
for  an  84  weight-percent  solids  loaded  PBAN  propellant,  A  maximum  tensile 
stress  criterion  appears  consistent  with  the  limited  data  for  the  tensile 
octant;  an  internal  friction  criterion  appears  gcncrically  applicable  in  the 
octants  where  tensions  and  compression.s  are  combined.  Superposed  hydro¬ 
static  p-'essurc  on  uniaxial  tension  caused  only  small  changes  in  strength  and 
elongation. 

Section  3.3  describes  the  results  of  the  limited  studies  and  experiments  that 
have  directed  attention  to  relationships  betv/ecn  crack  propagation  velocity 
and  propellant  physical  characteristics.  Good  agreement  v/as  found  between 
measured  and  predicted  crack  extensions. 


3.2  FAILURE  SURFACE  STUDIES 


Stress  and  strain  analyses  are  two  fundamental  requirements  for  determining 
a  rocket  motor's  structural  capabilities.  The  linear  theory  of  viscoelasticity 
has  proviced  a  useful  (though  not  always  precise)  guide  to  the  stresses  and 
strains  in  rocket  motors.  Prediction  of  failure  in  solid  propellant  grains, 
however,  has  been  dependent  upon  a  much  less  reliable  source.  In  most 
instances,  fa.lure  analyses  lean  heavily  on  uniaxial  tensile  test  data.  Since 
a  solid  propellant  grain  is  subjected  to  a  spectrum  of  three-dimensional 
stress  states,  the  applicability  of  uniaxial  tensile  data  to  the  prediction  of 
motor  failure  is  questionable.  An  understanding  of  how  failure  properties 
depend  upon  the  state  of  stress  or  stress  axiality  is  required  for  the  proper 
conduct  of  a  failure  analysis. 
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One  descriptive  correlation  of  failure  dependence  on  stress  state  is  the  so- 
called  failure  surface,  defined  in  a  three-dimensional  space  by  the  values 
of  the  principal  stresses  at  failure.  Th’  failure  surface  and  its  general 
topology  have  been  discus  sed  elsewhere  (Rtf.-.i),  A  major  difficulty  associ¬ 
ated  with  this  type  of  representation  is  the  fact  that  the  failure  properties  of 
solid  propellants  depend  upon  their  loading  history  as  well  a-s  stress  axiality. 
Since  a  given  stress  state  at  failure,  in  principle,  may  be  reached  by  an 
infinite  variety  of  loading  histories,  a  failure  surface  exists  only  for  specific 
loading  histories.  Sharma  (Ref.  a  ),  for  example,  describes  traces  of  the 
failure  surface  for  those  loading  histories  in  which  the  maximum  principal 
strain  rates  are  identical. 

Experimental  results  on  the  stress  state  dependence  of  filled  rubber  failure 
properties  have  been  reported  previously  by  Sharma,  Ref.  S  )  Jones, 

(Ref.  I  )  and  Majerus,  (Ref.  3  ).  In  Sharina's  paper  (Ref.  6  ),  attention 
was  focused  upon  delineating  the  failure  surface  for  a  viscoelastic  material, 
which,  however,  v/as  not  representative  of  most  propellants.  Jones  and 
Majerus  made  no  specific  attempt  to  establish  the  failure  surface. 

The  follov,'ing  considers  the  failure  surface  for  a  PBAM  propellant  with  a 
solids  loading  of  between  69  and  80  percent  by  volume.  The  mechanics  of 
the  test  and  the  reduction  of  data  to  a  reasonable  presentation  are  empha¬ 
sized.  Representative  data  for  one  strain  rate  and  one  temperature  are  pre¬ 
sented  for  a  number  of  multiaxial  stress  states.  While  a  number  of  questions 
remain,  particularly  with  respect  to  proper  stress  analysis  of  the  tests,  the 
results  of  this  investigation  are  encouraging.  A  fir.st  approximation  to  the 
failure  surface  in  stress  space  has  been  obtained,  and  its  variation  with  pro¬ 
pellant  solids  loading  has  been  partially  explored. 


3.2.1  Test  Analysis 


While  linear  viscoe'astic  stress  analysis  has  provided  a  significant  guide  to 
propellant  stress  analysis,  it  po.sstsses  limitations  which  are  not  fully 
understood,  e.g.  ,  large  deformation  effects  and  the  differing  responses  to 
compressive  and  tensile  loads.  Although  it  is  necessary  to  consider  these 
effects  if  an  analysis  of  failure  data  is  to  be  precise,  they  cannot  be  accounted 
for  at  the  present  time  in  a  refined  manner.  Indeed,  their  incorporation  in 
test  data  analysis  depends  as  much  upon  engineering  extrapolation  as  upon 
mathematical  rigor.  As  a  consequence,  the  following  work  has  to  be  viewed 
more  as  a  step  in  the  proper  direction  than  as  an  ultimate  result. 

For  the  purpose  of  the  present  discussion,  we  shall  distinguish  between 
effects  of  large  deformations  upon  the  stress  field  and  the  effect  of  nonlinear 
material  behavior.  It  can  be  shown  (Ref.  6  ,  7  )  that  the  change  in  stress 
values  due  to  consideration  of  large  deformation  rarely  exceeds  the  order  of 
the  maximum  strain  value.  While  deformations  at  failure  in  dominantly  com¬ 
pressive  stress  fields  may  exceed  25  percent,  deformations  in  the  tests 


3-2 


LOCKHEED  PROPULSION  COMP 


repor ied  here  are  generally  of  this  order  or  less.  The  stress  analyseE  of 
the  tests,  therefore,  are  bac-ed  on  the  usual  assumption  of  infinites Imal 
deformations.  The  effects  of  bilinear  cornpres sion-tension  mod;}U  are  dis¬ 
cussed  in  the  lollov/ing  paragraphs  appropriate  to  tl  e  particular  tests  in  v/aich 
they  are  an  important  factor. 

Seven  testa  reprt fienting  e  multiplicity  of  locations  in  three -dime nsional 
etresfj  space  have  been  in  compiling  the  data  presented  herein.  Test 

specimen  geometries  are  shown  in  Figure  3-1.  The  tests  are  listed  below, 
and  their  principiil  stress  values  are  shov/n  for  the  observed  sites  of  failure 
initiation  in  terms  of  applied  boundary  stresses  or  loads.  The  calculation 
of  these  stresses  is  based  on  the  assumption.?  of  material  linearity,  homo¬ 
geneity,  continuity,  and  infinitesimal  deformation.  Arabic  subscripts 
indicate  principal  stresses,  the  "1"  direction  identifies  the  principal  tensile 
stress  in  each  test. 

Strip  Biaxial  Tension 


(r3  =0 

where  co  is  the  average  applied  stress  for  a  strip  of  length  at  least  5 
times  its  height  (Figure  3-1). 

Uniaxial  Tension 


<r,  =  (To 

o’i  =  ffj  =  0 

where  is  the  applied  tension. 
Uniaxial  Compression 

0-1  = 

(Tj  =  <r,  =  0 

where  cr^  is  the  applied  compression. 
Strip  Biaxial  Compression 

0-1  =  -<ro 

<^2  =  <^1/2 


(2) 


(3) 


(4) 
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where  ijs  the  applied  cornprcftciCii  for  &  etrip  cf  length  St  Ic&st  S 
time®  Its  hsi'ht, 

V-' 


D  i.  T;  m  r.  t  r  s  1  Co ’t-  n  r  o  *1 i  o  n 

F 

e  1  " 

•  V  <  <  Ci  i 

^2  =  -BcTj  (5) 

e,  =  0 
where 

F  =  applied  compression 
d  =  specimen  diameter 
t  =  specimen  thickness. 

Shear 


ff,  =  T 

ffj  =  -To  (6) 

e,  =  0 

where  Is  the  average  applied  shearing  stress  in  a  laterally  restrained 
or  "simple"  shear  geometry  (Figure  3-1). 

Uniaxial  Tension  V/ith  Superposed  Hydrostatic  Pressure 


where  -tr  ^  Is  the  applied  hydrostatic  gas  pressure. 


(7) 


The  shear  test  and  the  diametral  compression  test  require  additional  con¬ 
siderations.  Previously,  it  was  determined  that  shear  strains  produce 
significant  stresses  normal  to  the  applied  shear  (Ref.  9  )  if  the  outer  sup¬ 
ports  are  restrained  from  la^’eral  motion.  Specifically,  it  was  found  that 
the  normal  stress  e  y  across  the  shear  is  proportional  to  the  applied  shear 
stress  and  is,  approximately, 

(T  y  S  0. 3  (8) 

where  is  the  shear  force  divided  by  the  total  area. 
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Herrmann  (Ref.  lo)  has  considered  the  effect  of  bilinear  material  properties 
on  the  simple  shear  test  and  shows  that  the  normal  stress  ^  along  the 
largest  dimension  is  equal  to  cr  y  For  a  ratio  of  tensile -to-compre s sive 
modulus  of  1/2,  which  is  characteristic  of  propellants  with  the  solids  lond- 
ing«  investigated  here,  the  stresses  cr  ^  and  <t  y  have  been  found  to  be  approKi- 
mately  1/3  the  shear  stress.  Therefore,  the  principal  stresses  are 


OTj  =  0.7  Tj, 


(9) 


whereas  linear  theory  yielded  ctj  =  -o-^  =  (Equation  6). 


A  similarly  detailed  analysis  for  the  diametral  compression  test  has  not  been 
performed,  primarily  because  the  stress  distribution  varies  through  the 
body,  which  significantly  complicates  the  analysis.  However,  an  estimate 
of  the  correction  can  be  obtained  if  one  considers  the  different  modulus 
values  in  tension  and  compression.  Inasmuch  as  the  compressive  modulus 
is  about  twice  as  large  as  the  tensile  modulus,  one  could  assvime  that  half 
of  the  compressive  force  is  absorbed  in  the  disk  without  producing  trans¬ 
verse  stresses  and  that  the  remainder  of  the  force  gives  rise  to  the  tensile 
stress  (Tj.  Instead  of  the  31  ratio  between  the  principal  compression  and 
tension  stresses  found  in  linear  analysis  (Equation  5),  a  ratio  nearer  to  6:1 
appears  to  be  more  realistic. 


It  is  clear  that  nonlinear  material  behavior  can  influence  significantly  the 
evaluation  of  the  test  results.  While  the  precision  of  the  analysis  is  question¬ 
able,  it  is  qualitatively  correct  and  reduces  the  effects  due  to  material  non¬ 
linearity. 


3.2.2  Experimental  Results 


Tests  were  performed  on  an  Instron  testing  machine.  The  testing  speed  was 
adjusted  so  that  the  principal  tensile  strain  rate  Cj  was  equal  to  0.07  in. /in./ 
min  when  based  on  the  specimen  dimensions  before  deformation.  The  occur¬ 
rence  of  failure  was  defined  by  observed  cracking  or  by  the  maximum  in  the 
load -displacement  trace,  whichever  occurred  first. 

The  test  points  shown  here  are  representative  of  a  more  complete  set  of 
failure  properties  given  in  Reference  1.  For  the  present  purpose 

of  examining  the  effect  of  stress  axiality  on  failure,  it  suffices  to  consider 
data  derived  from  one  strain  rate  and  at  one  temperature.  Figures  3-2  and 
3-3  show  failure  data  for  three  solid  loadings  of  a  PBAN  binder  at  10®F  and 
7  percent/min  strain  rate.  The  shear  and  diameiral  compression  test  data 
in  Figure  3-3  w'cre  reduced  according  to  the  bilinear  analysis. 


This  ratio  depends  also  upon  loading,  but  such  detailed  consideration  has 
been  omitted  here. 
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3.2.3  Discussion 


The  modes  of  failure  observed  v.'ere  normal  tensile  cracking  in  the  strip  and 
uniaxial  tensile  tests  and  shear  sliding  for  the  remaining  tests.  Shear  fail¬ 
ure  surfaces  were  oriented  approximately  30”F  off  tVie  direction  of  the  rnaxi- 
rnum  normal  compression  etress.  An  obvious  increase  in  the  shear  stresses 
accompanies  increases  in  normal  compressions.  This  is  believed  to  be  a 
consequence  of  intergranular  friction  under  compression  and,  of  course,  is 
not  significant  under  tensions  where  the  granules  tend  to  be  separated  by  the 
deformations.  In  this  respect,  the  general  failure  behavior  of  the  propellant 
is  similar  to  that  of  a  water  -  saturated  sandy  clay  (Ref.  il).  Such  materials 
may  be  characterized  by  a  single  value  of  tensile  strength  and  a  shear  strength 
that  varies  in  proportion  to  the  normal  stress  across  the  plane  of  potential 
shear  fracture.  This  kind  of  behavior  is  represented  conveniently  in  the  form 
of  a  Mohr  envelope  (e.  g.  ,  Ref.  12). 

Uniaxial  data  obtained  with  tension  superposed  on  a  hydrostatic  (gas)  pre.s- 
sure  for  the  70  percent  by  volume  solids  loaded  formulation  showed  that  the 
amount  of  tension  required  to  cause  failure  increased  only  slightly  for  pres¬ 
sures  up  to  1000  psi  (Figure  3-2).  The  indicated  tensile  strength  is  nominally 
con  tant  and  independent  of  the  pressure,  any  increment  of  the  latter  being 
cariied  by  neutral  pressures  in  the  rubber  matrix  which  fills  the  voids 
between  particles.  These  neutral  pressures  in  the  rubber  naturally  do  not 
contribute  to  an  increase  of  the  intergranular  frictions  or  binding  forces 
between  the  granule?:  and  the  matrix. 

The  small  increa.ses  noted  in  strength  with  pressure  can  be  attributed  to 
increases  in  intergranular  friction,  interphase  friction  or  increases  in  the 
strength  of  the  binder  alone  under  pressure.  The  size  of  the  effect  is  such 
that  any  of  the  three  phenomena  could  rationalize  the  strength  changes;  as 
the  solids  volumetric  loading  is  increased,  the  interparticle  reactions  can 
be  expected  to  dominate  in  importance. 

Since  each  of  the  stress  states  examined  in  these  tests  (excepting  pressurized 
tensile  tests)  lies  on  one  of  the  coordinate  planes*  in  the  principal  stress 
space,  an  experimental  picture  of  the  complete  failure  surface  cannot  be 
constructed.  However,  in  view  of  the  need  of  failure  analysis  in  solid  pro¬ 
pellant  rocket  design,  it  becomes  interesting  to  sketch  a  possible  extrapolation 
of  the  biaxial  data  to  the  triaxial  stress  states.  Tests  to  examine  the  failure 
surface  in  these  octants  are  planned  for  the  future.  As  a  practical  example 
of  such  need  for  extrapolation,  consider  the  bond  stresses  in  an  unpressur¬ 
ized,  axially  accelerated  motor;  in  this  case,  the  principal  stresses  can  lie 
in  the  +  +  -  or  +--  octants,  substantially  away  from  any  of  the  coordinate 
planes.  Extrapolation  of  the  data  of  Figure  3-3  to  construct  a  reasonable 
surface  may  be  done  in  various  ways.  Figure  3-4  illustrates  a  construction 
for  the  part  of  the  sur''  <-  above  the  o-j,  plane,  based  on  use  of  the  Mohr 
envelope . 

*  It  should  be  recalled  '-’at  we  are  considering  engineerii;g  stresfc  ’s  only.  Due 
to  large  deforrnaliun?i .  Uie  stress  state  may  be  charged  such  that  the  test 
points  do  not  fall  exac.ly  on  the  coordinate  planes. 
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The  data  obtainrd  in  this  study  do  not  fit  any  of  the  classical  arialytic  failure 
theories  with  convincing  accuracy.  At  beet,  the  limited  data  on  the  planes 
in  the  +++  octant  suggest  a  rnaximvim  tensile  stress  criterion,  and  tb,c  elicar 
and  ccsnaprcEsioi'i  data  extrapolsded  into  the  +--  octant  arc  not  in.  gross  dis¬ 
agreement  v.'ith  a  linear  internal  frictio.i  criterion.  However,  agreement 
Vi;ith  an  analytically  formula.ted  criterion  ia  valuable  but  not  critical  for 
design  purposes  and  this  initial,  pra.grnetic  ntapping  of  the  surfaxe  produces 
in  principle  the  desired  engineering  inforination.  With  connective  theory 
relating  the  mechanics  cf  failure  to  the  shape  of  the  failure  surface,  a  con¬ 
fident  reduction  in  engineering  detail  required  for  the  mapping,  in  addition  to 
an  insight  into  material  st-ength  optimization,  would  result.  This  is  an 
ultimate  objective  of  the  work  in  progresB. 

Experimental  tests  of  PBAA,  carboxy-terminated  and  Nitroplastisol  propel¬ 
lants  in  compression,  shear  and  tension  have  indicated  that  the  general  failure 
behavior  of  these  propellants  is  not  widely  at  variance  with  that  depicted  in 
figures  3-3  and  3-4.  Of  course,  there  are  differences  in  stress  magnitudes 
at  failure  and  in  the  envelopes  for  different  strain  rates  and  temperatures. 

It  is  noteworthy  that  the  shape  of  the  failure  surface,  as  well  as  its  spatial 
location,  can  change  significantly  with  temperature  for  a  single  formulation, 
of  the  propellant  (e.  g,  ,  Ref.  1,  ll). 

It  should  be  kept  in  mind  that  the  above  data  and  their  interpolation  lack  the 
support  of  a  precise  stress  analysis.  Not  withstanditig  this  deficiency,  the 
information  as  presented  exemplifies  a  better  engineering  tool  for  rocket 
motor  evaluation  than  previously  has  been  available.  The  obvious  limitations 
of  the  data  emphasize  the  need  for  additional  experimental  points  on  the  failure 
surface  off  the  coordinate  planes  and  practical  stress  analysis  methods  that 
account  for  intrinsic  and  deformational  nonlinearities  in  propellant  behavior. 
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3.3  FRACTURE  PROPAGATION  STUDIiCS 


Fracture  propagation  in  viscoelastic  materials  is  regulated  by  their  strength, 
the  driving  stresses  and  the  energy  absorbed  in  the  fracture  process.  Each 
of  these  parameters  is  time  and  temperature  dependent.  This  section 
describes  the  results  of  a  study  of  crack  propagation  in  biaxially  stressed 
plates  of  PBAA  propellant. 

Current  theories  of  fracture  in  polymers  gcne’*ally  agree  that  crack  e.xten- 
sions  occur  in  discrete  jumps  of  a  distance  6  which  occur  in  a  time  t6  .  6  is 

a  fixed  dimension  characteristic  of  the  material.  It  is  proposed  herein  that 
t6  is  a  particular  function  related  to  the  measured  failure  properties  of  the 
propellant  and  the  driving  stress.  The  time  dependence  of  (he  fvanction  at 
various  temperatures  is  derived  and  compared  with  experimental  data  for 
crack  extension  and  crack  velocity  under  conditions  of  time -varying  driving 
stress.  Excellent  agreement  between  measured  and  predicted  crack  exten¬ 
sions  was  observed. 


3.3.1  Analysis 


The  first  criterion  for  the  elastic  fracture  of  rubbery  materials  resulted 
from  an  extension  of  the  Griffith  brittle  fracture  theory  (Ref.  13)  by 
Rivlin  and  Thomas  (F..ef.  14).  Greensmith  and  Thomas  (Refs.  15,  16,  17,  18) 
observed  that,  under  equilibrium  conditions,  the  energy  absorbed  by  fracture 
could  be  approximated  by 

T  =  Wjd  (1) 

where  Wj  is  the  critical  strain  energy  density  for  failure  in  an  initially 
uncracked  specimen  and  d  is  twico  the  radius  of  curvature  of  the  crack  tip, 
Williams  (Refs.  19,  ?0)  applied  the  concept  of  the  Griffith  criterion  to  predic¬ 
tion  of  crack  initiation  and  propagation  in  rubbery  materials  along  similar 
lines  of  inquiry. 


Williams'  rnodrl  of  the  crack  postulated  that  the  crack  extended  in  discrete 
jumps  of  distance  6  .  In  an  infinite  sheet  stressed  biaxially  in  tension,  the 
stress  at  the  tip  of  a  crack  of  initial  length  2b  in  the  sheet  is  given  by 


7J- 


b  +  m6 

21 - 


(2) 


where  is  the  stress  normal  to  the  crack's  axis  at  a  distance  and  m  is  the 
number  of  jumps  by  which  the  crack  has  extended.  The  crack  tip  radius  of 
curvature  d/2  is  identified  in  Equation  (2)  as  6.  According  to  this  relation¬ 
ship,  the  stres.s  increases  with  crack  length  b  +  m6and  the  crack  velocity  may 
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be  expected  to  incrcncc  v/ith  its  extension.  Applying  linc.-?.r  tlieory  to  t’<,e  cr.-;.ck 

in  a  viecoclastic  medium,  V/illiarna  (P.ef.  20)  found  crack 

velocity  and  extension  with  time.  These  expressions  were  based  on  the  assurnp 

tion  of  a  maximum  strain  rupture  criterion. 

Recent  work  at  LPC  (Refs.  21,  1,  22)  has  eb'wn  that  gross  failure  in  prc]:*?:!- 
lant  conforms  to  a  maximum  tensile  stress  iailure  criterion.  Application  of 
the  stress  criterion  to  the  crack  extension  problem  in,  tlscrcfore,  apu'iropriate. 

Only  a  special  case  of  crack  extension  will  be  treated  bare:  namely,  when 
the  crack  tip  stresses  arc  not  a  function  of  crack  length.  Equation  (2)  above 
describes  the  stresses  at  the  crack  tip  when  the  crack  is  small  with  respect 
to  the  size  of  the  plate  and  far  removed  from  its  boundaries.  When  the  crack 
is  long  with  respect  to  the  sheet  width  (perpendicular  to  the  crack),  the  crack 
tip  stresses  do  not  change  materially  with  extension.  In  practice  this  condi¬ 
tion  is  achieved  by  clamping  a  long  plate  between  rigid  boundaries  and  slitting 
it  a  distance  about  equal  to  its  height  (or  gage  length).  The  slitted  specimen 
is  then  extended  normal  to  the  slit,  which  can  be  considered  to  be  an  initial 
static  crack.  The  critical  dimensions  of  the  crack  inflviencing  its  propagation 
under  these  boundary  conditions  are  its  tip  radius  and  the  jump  distance  6. 

The  crack  extension  model  to  be  considered  here  follows  that  proposed  by 
Rivlin,  Thom.as,  Williams,  et  al,  in  that  extension  of  the  crack  occurs  in  dis  ¬ 
crete  jumps  of  distance  6.  When  the  extension  4  occurs  in  a.  time  t6  ,  the 
crack  velocity  S  is  given  by  the  identity 


The  variation  of  S  with  material  propertiec  and  crack  driving  stress  is  .speci¬ 
fied  by  the  variation  of  t6  with  these  properties.  6  is  a  constant  characteristic 
dimension  for  the  material  thought  to  be  related  to  its  microstructure. 

Extending  consideration  of  Equation  (3)  for  the  viscoelastic  case,  we  postulate 
that  S  (t)  is  proportional  to  the  tirre  to  failure  observed  in  tests  of  the  material 
and  inversely  proportional  to  the  driving  stress.  The  particular  test  mode  of 
interest  herein  is  constant  strain  rate.  When  measured  at  constant  strain 
rate,  the  rupture  stress  in  propellant  varies  with  the  time  to  rupture  of  the 
sample  as  shown  in  P'igure  3-5.  The  observed  relationship  may  be  approxi¬ 
mated  analytically  by 

0^(t)  ft.  B(l/t)"^  (4) 

over  wide  ranges  of  time.  Assuming  gross  rupture  behavior  varies  identi¬ 
cally  with  the  microscopic  fracture,  the  variation  of  tfl(t)  with  time  may  be 
expressed  as 

t6{t)  (l/t)”^  (5) 
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where  kj  is  an  urdeterrained  constant.  The  driving  strccs  in  the  sheet  under 
constant  strairi  rate  cCinditions  maybe  apprcsiirurited  by 


it) 


fcj 


Af  .1-n 


1-n 


(6) 


when  the  relaxation  mcdvdiis  is  linear  with  tuT.e  in  log-log  ecordinates.  n  is 
the  slope  of  t’le  log  inodvd-us  curve,  A  is  its  intercept  at  t  -■  I,  and  s  t!:is 
©train  rate.  The  rclajcnncn  modulus  for  the  propellant  used  in  these  experi¬ 
ments  is  shov;n  in  Figure  3-6. 


Postulating  that 


t6(t)  w 


k 


B(l/t)'"" 


it  follows  that 


.  .  .  6A* 

s(t)^  ftt^T 


1  -(m+n) 


and,  integrating 


(7) 


(8) 


kj  6  A< 


^2-(m+n) 


(9) 


S(t)  plotted  versus  t  should  yield  a  linear  relationship  for  the  crack 

propagation  model  deecr.bcd  by  Lquatioa  (9).  The  crack  enicnsion  with  time 
is  an  experimental  observable. 


According  to  the  work  of  Williams  (Ref.. 20), the  dimension  may  be  estimated 
as  the  critical  radius  of  the  crack  tip.  The  measured  critical  crack  tip  dia¬ 
meter  is  approximately  0.02  inch  and,  hence,  6  *a  0.01  inch.  The  value  of  6 
for  gum  rubbers  was  estimated  to  be  of  this  order  by  Thomas  (Ref.  17).  The 
constants  kj  and  c  are,  therefore,  the  only  undetermined  constants. 


3.3.2  Experimental  Results 


Fracture  extension  experiments  were  made  using  biaxial  strip  test  specimens 
of  PBAA  propellant  (LPC-543B,  Batch  4404).  The  specimens  were  machined 
plates  (10  X  2.25  x  0.13  inches)  within  situ  bonded  wooden  rest.aints  along 
the  long  edges.  Fracture  specimens  were  slit  v/ith  a  razor  blade  along  their 
midlines  a  distance  of  2.4  inches  from  one  cdg<  .  The  specimens  were  tested 
in  an  Instron  tester  at  a  ct.  sshead  speed  of  0.1  inch/minute  to  obtain  a  strain 
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igure  3-6  Master  Sti'ess  Relaxation  Modulus  and  Shift  Factor, 
Propellant  L.PC-543B  (}3atch  4404) 
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rate  of  0,0444  mln“*  in  the  body  of  plate.  The  ultimate  property  data  gho%vn 
in  Figure  3-5  were  obtained  in  tests  of  identical  but  unfcUt  specivaons  @.t  the 
same  strain  rate.  The  stress  relaxation  modulus  data  shown  in  Figure  3-6 
were  obtained  in  tests  of  uniaxial  wooden  tab-end  tensile  Epc.cin;''n5:.  t*Hsted 
in  accordance  with  Iht  ICRPG  test  procedure  (CPIA  No.  El,  PFSS).  Tise 
data  of  Figures  3-5  and  3-6  are  prcEcnted  in  reduced  variable  form,  \ojing 
the  shift  factor  shown  on  Figure  3-6. 

Crack  extensions  were  measured  at  0.1  inch  intervals  by  visual  observation 
of  the  crack  through  a  magnifying  glass.  Typical  data  are  shown  in  Figure  3-7. 
Measured  crack  extensions  are  show'n  plotted  versus  in  Figure  3-8. 

The  exponents  n  and  rn  were  measured  from  shifted  data.  The  rdtopcs  were 
taken  in  the  time  region  within  which  the  experiment  was  done  and  the  inter¬ 
cepts  at  various  temperatures  determined  graphically. 


3.3.3  Discussion 


The  data  of  Figure  3-8  show  excellent  agreement  with  the  postulated  time 
variation  of  crack  extension.  Taking  the  constant  of  integration  in  Equation  (9) 
as 

,  eA  i 

~  T<.jB(  1  -n)(E-n-m) 

where  t|  is  the  time  of  crack  initiation  and  t  -  0.01,  the  undetermined  constant 
kj  is  found  to  be  165  for  the  data  presented.  Crack  velocities  predicted  by 
Equation  (8)  are  of  the  order  of  1.0  to  2.5  inches  per  minute,  depending  on 
crack  length  and  test  temperature,  in  excellent  agreement  with  the  experi¬ 
mental  observations. 

The  agreement  between  the  crack  velocity  and  extension  relationships, 
Equations  (8)  and  (9),  apply  only  to  the  constant  strain  rate  test  mode  for 
initially  quite  long  cracks.  Further  analysis  is  required  for  this  test  mode 
for  cracks  initially  small  with  respect  to  plate  dimensions.  Small  cracks 
can  be  expected  to  accelerate  rather  rapidly  with  time,  because  of  the  high 
stress  concentrations  at  their  tips  (cf.  Equation  2). 

The  failure  properties  of  propellant  (i.  e.  ,  those  shown  in  Figure  3-5)  do 
not  shift  linearly  with  varying  strain  rate.  Extension  of  S(t)  and  S(t)  relation¬ 
ships  to  consideration  of  crack  extension  at  varying  strain  rate  requires, 
therefore,  measurement  of  the  failure  properties  at  various  strain  rates. 

Crack  propagation  under  other  coi,ditiona  of  time  varying  stress  is  a  matter 
of  both  gt  neral  and  specific  interest.  Tiie  observed  slope.s  of  rupture  stress 
versus  time  curves  for  constant  stress  and  constant  strain  tests  of  the  pro¬ 
pellant  used  in  these  tests  is  approximately  the  s.ame  as  the  slope  of  the 
curve  shown  in  ~igure  3-5.  a  he  ar.alysis  for  Si't  t  appropriate  to  the  constant 
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Figure  3«7  Crack  Extension  versus  Time 
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BtrcBs  test  condition  predicts  that  the  cracks  propagate  at  a  nearly  uniform 
velocity.  Under  constant  strain  conditions,  a  crack  started  by  the  initial 
straining  should  decelerate  and  stop.  Both  predictions  are  in  accord  v/ith 
data  obtained  in  single  tests.  Further  analysis  and  experimentation  defining 
such  behavior  is  rccorr.n'.cndcd. 

The  agreen-ier.t  between  the  analytical  predictions  and  experimental  observa¬ 
tions  Is  conGidered  excellent.  Extended  efforts  v.dth  greater  attention  to 
rigorous  theoretical  development  of  tiie  analysis  is  recommended.  Future 
engineering  applications  of  fracture  mechanics  analysis  include,  for  example, 
rocket  grain  defect  analysis  and  prediction  of  threshold  conditions  for  catas¬ 
trophic  grain  cracking.  Extension  of  fracture  mechanics  analysis  to  the 
energetics  of  fracture  extension  is  also  a  potentially  fruitful  method  for 
studying  the  relationships  among  binder  molecular  structure,  interphaoe 
adhesions,  solids  particulate  properties,  and  intrinsic  propellant  strength. 
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4.1  INTRODUCTION 


Current  areas  of  interest  arc  as  follo-.vs: 

0  Determination  of  viscoelastic  and  thermoelastic  properties  of 
propellant  and  associated  materials,  using  cyclic  techniques 

0  Investigation  of  propellant  mechanical  and  thermal  response  to 
sustained  cyclic  loading 

0  Determination  of  vibration  failure  modes  in  solid  propellantE 

and  correlation  of  propellant  microstructure  with  vibration  failure 
resistance  ^ 

Previous  w.ork  at  LPC  (Ref.  21  )  has  demonstrated  the  equivalence  of  transient 
and  dynamic  physical  property  characterization  techniques  for  practical  pro¬ 
pellants.  The  problems  introduced  by  the  nonlinear  viscoelastic  behavior  of 
highly  loaded  propellant  systems  (c.g.  ,  strain-induced  nonlinearity  and 
thermal  degradation)  v.'cre  examined  Iny  means  of  large  deformation  dynamic 
tests.  In  addition,  viscoelastic  heating  effects  were  studied;  and  the  possi¬ 
bility  of  regenerative  he.ating,  leading  to  thermal  and  raachar.ical  instability 
in  solid  propellant  systems  subjected  to  sustained  cyclic  loading,  v/as  deduced 
on  theoretical  grounds.  The  current  work  is  directed  at  the  extension  and 
amplification  of  these  results. 

Experimental  investigations  are  discussed  in  which  propellant  dynamic  shear 
and  bulk  properties  were  measured  with  the  small  deformation  piezoelectric 
devices  previously  described  (Ref.  1  ).  Shear  moduli  of  PBAA  and  PBAN 
propellant  and  cured  binder  are  presented  and  compared  with  linear  thermal 
expansion  measurements  of  glass  transition  temperature  obtained  for  similar 
propellants.  Analysis  techniques  and  calibration  results  are  discussed  for 
the  dynamic  bulk  modulus  device. 

An  experimental  investigation  of  propellant  thermomechanical  response  to 
sustained  cyclic  inertial  loading  was  completed.  Experimental  results  are 
presented,  noting  that  they  are  in  qualitative  agreement  with  theoretical  pre¬ 
dictions.  In  particular,  jump  instability  caused  by  regeneral.ve  thermal 
and  mechanical  coupling  was  obtained  experimentally.  Also  completed  were 
experimental  investigations  of  the  transient  thermoviscoelastic  response 
under  constant  cyclic  strain  amplitude  (fixed  output)  loading  and  with  inertial 
loading.  Experimental  results  are  given  including  an  inve stigation  of 
cyclic  degradation  effects. 
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deformatlcn  a.re  measured  and  arc  used,  together  Vvitli  the  rAcasured  speci¬ 
men  dimensions,  to  calculate  dynamic  viscoelastic  propcrtieB  of  the  Kpccimen. 
The  apparatus  is  presently  being  used  to  viocoelasticElly  characterise  a 
wide  variety  of  propellants,  binders  snd  asBOciated  rocket  motor  materials. 


Dynamic  shear  modvili  of  PBAA  and  PBAN  propellants  and  cured  binders 
have  been  obtained  over  a  wide  range  of  temperature  and  at  frequencies  from 
20  to  iOOO  cycles  par  second.  Stor8m.e  and  loss  mcdvrli  (C  and  G",  respec¬ 
tively)  are  ehov/n  plotted  versus  temperature  at  several  frequencies  in 
Figures  4-1  and  4-2  for  both  the  S4  percent  solids  loaded  propellants  and 
the  cured  binders.  In  addition,  the  "glass  transition  temperature  region" 
as  measured  by  linear  thermal  expansion  coefficient  techniques  on  similar 
propellant  formulations  is  shown  on  each  graph.  Correlation  of  this  tem¬ 
perature  v/ith  the  temperature  region  v/here  the  propellant  dynamic  properties 
become  relatively  independent  of  temperature  and  frequency  was  observed  for 
Nitroplactisol  propellant  and  binder  in  addition  to  the  rubber -ba»..e  propellants 
shown. 

The  data  show  the  generally  anticipated  effect  of  the  addition  of  solid  filler 
to  an  elastomer.  The  viscoelastic  transition  region  is  significantly  broadened 
and  the  moduli  are  Increased  particularly  at  higher  temperatures.  Charac¬ 
teristic  regions,  such  as  the  loss  moduli  peaks  and  the  low  temperature  limit 
where  the  response  becomes  frequency  independent,  arc  ob.served  to  shift 
slightly  to  higher  temperatures  with  the  addition  of  filler.  This  effect  also 
occurs  for  the  Nitroplastisol  propellant  and  binder  v/here  the  magnitude  of 
the  shift  appears  to  be  much  more  significant.  The  anticipated  shift  of  the 
response  curve  to  higher  temperaturee  as  the  excitation  frequency  is  increased 
is  clearly  evident  for  both  sets  of  data.  The  generally  lower  characteristic 
temperatures  for  the  PBAA  as  compared  with  the  PBAN  materials  correlates 
well  with  the  better  low  temperature  physical  property  capability  observed 
for  PBAA  propellants. 
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Figure  4-1  Dynamic  Shear  ModviluE  versus  Tempernture,  PBAA 
Propellant  and  Cured  Binder,  84^  To"ai  Solids  Loadinj 
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Figure  4-2  D\  namir  Shear  ?vIodul  le  versus  1  e  tnpe  ratur  e,  PBAN 

Piopellarit  nrd  Cured  Binder,  84“^  I'otal  Solids  Loading 
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4.2,2  Dynimic  Bulk  Modulus 


An  apparatus,  similar  to  the  dynamic  shear  modulus  device  (described  in 
the  last  section)  but  designed  to  measure  t’  e  dynamic  bulk  moduli  of  materials, 
is  currently  being  evaluated.  The  apparaius  contains  a  pres.surc  cavity 
within  which  the  specimen  material  is  placed.  The  specimeii  ' hydrostatic¬ 
ally  pre  ssurii’.ed  and  a  small  sinusoidal  prcosu’c  is  superpo-sci  on  the 
hydrostatic  component  by  sinusoidally  varying  the  volume  of  the  cavity  by 
means  of  a  piezoelectric  driver  and  piston  arrangement  (Ref.  1  ).  The 
resulting  sinusoidal  pressure  is  measured  and  compared  with  the  input  vol¬ 
ume  change  to  compute  bulk  modulus  properties  of  the  confined  material. 

The  device  is  presently  being  calibrated  with  standard  fluids  for  eventual 
measurement  of  propellant  properties.  The  calibration  techniqtie  is  as 
follow's.  The  bulk  modulus  (K)  is  defined  as 

K  -  (1) 

where  V  is  the  nominal  cavity  volume,  A  V  is  the  dynamic  volume  change  and 
AP  is  the  resulting  dynamic  cavity  pressure.  Of  major  concern  in  the  cali¬ 
bration  process,  particularly  for  high  bulk  modulus  materials,  is  the  effective 
bulk  modulus  or  compliance  of  the  constraining  cavity.  Clearly,  for  a  real 
cavity  with  finite  compliance,  the  input  volume  change  will  equal  the  sum  of 
the  sample  volume  change  and  the  cavity  volume  change  ^  that  i.s: 

AV--AV_^+^Vq  (2) 

Since  both  the  cavity  v.'alli;  and  the  constrained  £an->ple  experience  the  same 
pressure  variation  £iP  (assuming  that  measurerrient  frequencies  are  low  enough 
to  make  pressure  wave  propagation  effects  negligible)  w’c  can  define  the  sam¬ 
ple  bulk  modulus  K|  and  the  effective  cavity  bulk  modulus  as  follows: 

.  Kg  =  (3) 

Combining  these  definitions  with  the  volume  relationship.  Equation  (2), 
results  in 

=  B  (4) 

where  B^  and  B^  are  the  sample  and  cavity  compliances,  defined  as  recip¬ 
rocals  of  the  respective  bulk  moduli  and  B  is  the  total  measured  compliance, 
the  reciprocal  of  K  in  Equetion  (1).  The  output  dynamic  pressure  AP  and 
input  dynamic  "olume  change  are  lincarily  related  to  the  apparatus  output 
and  input  transducer  voltages  v^  and  as  follows: 

LP  C,  Vg  .  AV  - 
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Therefore,  for  coe.stant  inp\it  voltage,  vq,  and  constant  cavity  conipliance, 

Bq,  the  sample  compliance  is  linearly  i  elated  to  the  reciprocal  of  the 
measured  outpvit  voltage.  Thus  for  a  constant  cavity  compliance,  tlie  appa¬ 
ratus  can  l>c  calibrate<l  liy  defining  the  linear  relation.ship  of  Equation  (5) 
from  rnea  sur  emer  t  s  of  samples  v.ith  known  bulk  compliance. 

A  typical  calibration  curve  is  sho'wn  in  Figure  4-3  for  a  single  temperature 
and  hydrostatic  pressure.  In  general,  the  calib’-ation  curve  will  Vie  a  function 
of  both  temperature  and  hydrostatic  pressure,  the  temperature  dependence 
as.soc;iated  primarily  with  the  input  and  output  transducer  coefficients  and 
the  pressure  dependence  associated  primarily  with  a  reduction  of  mechanical 
slack  in  tlie  system  as  the  hydrostatic  pressure  is  increased. 
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4.  3 


LARGE  STRAIN  LYN 
HEATING  EFFEC'IS 


MIC  i'ROPERTlES  AND  VISCOELASTIC 


Earlier  work  at  I. PC  in  ‘.h-.s  area  con«^icted  of  analytical  studies  of  the  physical 
and  thermal  rosionse  ol  propellant  to  dynamic  loading  a;)d  experimental  meas¬ 
urements  of  the  riynarnic  rcstiorise  of  propellants  to  Lirge  cyclic  deformations 
(Refs.  2  1,  t  ).  These  investigations  showed  that  o.  large  temperature  rise 
could  be  obtained  an-,  ertain  types  of  lonlinearities  and  degradation  effects 
were  found  to  be  sig.  ’.icant  in  the  response  of  an  84  percent  total  solids-loaded 
propellant  at  strain  levels  of  'nly  a  few  percent. 

Present  work  is  exar.iining  the  validity  of  u.sing  linear  viscoelastic  theory  to 
predict  the  thermo-mechatiival  response  of  highly  filled  propellants.  The 
ultimate  aim  is  to  develop  engineering  methods  for  ihe  prediction  of  uyiiami>^ 
behavior  and  heating  effects  in  real  motor  situations. 

A-dditional  experimental  data  are  given  fioni  constant  strain  amplitude  and 
inertial  loading  experiment. s  on  propellant.  Regenerative  coupling  between 
generated  heat  and  viscoelastic  properties  was  obtained  under  inertial  load¬ 
ing  conditions  and  re.sulted  in  jump  -  ms  tab:  lity ,  a.s  predicted  theoretically. 

Further  experimental  data  are  presented  from  tests  in  which  transient, 
dynamic  the’-movl,sc.oela.stic  response,  and  degradat ion  effects  were  studied 
under  two  types  of  therrn.'il  boundary  conditions:  (a)  approximately  adiabatic 
.‘irui  (b),  unidirectional  heat  flow  to  a  constant  temperature  boundary. 

4.  3.  1  Experimental  Techniques 


Propellant  specimens  were  subjected  to  steady -sta*e  sinusoidal  shear  dis¬ 
placements,  as  .shown  in  Figures  4- -I  and  4-5  .  Two  loading  conditions 

were  investigated:  (l)  output  rigidly  connected  to  force  transducer  (fixed 
output),  and  (Z)  output  inernatiy  ioauud  by  large  mass. 

For  the  first  loading  condition,  the  input  strain  and  the  output  force  were 
measured  together  with  the  phase  relationship  betvveen  them.  From  these 
data,  the  dynamic  shear  storage  and  loss  moduli  were  then  obtained  as  pre¬ 
viously  described  (Ref.  1  ).  When  an  inertial  loading  test  was  made,  the 
irput  strain  amplitude  (Vj)  and  the  amplitude  of  the  mass  m)  were  meas¬ 
ured  and  also  the  phase  difference  (a)  between  them. 

In  both  test  modes,  the  specimen  is  driven  at  a  constant  sinusoidal  displace¬ 
ment  by  means  of  a  mechanical  cam-drive  arrangement.  This  is  shown 
pictorially  in  Figure  4-5  ,  arranged  for  an  inci  ti^il  leading  test.  The  milling 
machine,  shovm  at  the  right  of  the  photograph,  actuates  a  cam  and  cam -follower 
arrangement  to  drive  the  propellant  shear  specimen  sliown  in  the  center  of  the 
photograph.  Return  force  for  the  cam -follower  la  supplied  by  the  compre.ssion 
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spring  between  the  driver  shaft  guides.  Input  motion  is  measured  at  the 
sample  supports  by  a  potentiometer -type  linear  motion  transducer.  The 
center  support  of  the  double -lap  shear  specimen  is  attached  to  a  loading 
shaft  constrained  to  move  only  along  the  axis  of  shear  by  a  set  of  rectangular 
linear  motion  bearings.  The  inertial  mass  is  clamped  to  the  loading  shaft 
and  is  adjusted  by  addition  or  removal  of  weights.  The  sliding  friction  of 
the  arrangement  when  fully  loaded  is  only  a  few  ounces.  Output  motion  is 
measured  at  the  sample  support  by  another  potentiometer -type  linear  motion 
transducer . 

For  the  fixed  output  test,  the  inertial  loading  system  is  replaced  by  a  rigid 
force-measuring  cell  to  which  the  specimen  is  attached  (Ref.  2  1  ). 

The  excitation  frequency  is  changed  by  varying  th  i  cam  revolution  speed  in 
discrete  steps  from  1.3  to  45  cycles  per  second.  The  input  displacement 
amplitudes  are  controlled  by  interchanging  driver  cams. 

Specimens  were  prepared  by  bending  milled  propellant  bars  to  the  support 
members  with  epoxy  adhesive.  The  semicircular  fillets  shown  in  the  figures 
(4-4  and  4-5  )  were  machined  into  the  propellant  bars  to  prevent  specimen 

cracking  in  the  regions  of  high  stress  at  the  specimen  corners  .  All  tests 
were  performed  with  the  outer  supports  constrained  to  pre\  ent  lateral  strain 
and  standardize  the  test  geometry. 

Two  thermal  boundary  conditions  were  investigated.  The  first  was  an  app¬ 
roximation  to  unidirectional  heat  flow  laterally  outward  from  the  double -lap 
shear  center  support,  and  the  second  was  an  approximation  to  adiabatic 
conditions.  Figures  4-4  through  4-5  illustrate  the  unidirectional  heat  flow 
experimental  arrangement.  In  application,  all  exposed  propellant  surfaces 
were  thermally  insulated  with  a  soft  foam-rubber  material  so  that  the  major¬ 
ity  of  the  heat  generated  in  the  propellant  was  conducted  laterally  away  from 
the  low  thermal  conductivity  center  support  member  to  the  metal  side  supports 
which  were  maintained  at  constant  temperature  by  circulating  fluid.  An  app¬ 
roximation  to  adiabatic  conditions  v/as  achieved  by  replacing  the  conductive 
metal  outer  supports  with  wood,  a  relatively  low -conductivity  material,  and 
insulating  all  surfaces  with  foam-rubber  material.  The  temperature  distrib¬ 
utions  throughout  the  specimens  were  monitored  with  small  thermocouples 
imbedded  directly  in  the  propellant  or  support  material. 


4.3.2  Dynamic  Inertial  Loading  Tests 


An  interesting  prediction  obtained  from  the  analysis  of  the  inertial  loading 
problem  was  that  regenerative  coupling  between  the  temper ature -dependent 
properties  and  the  generated  heat  could  lead  to  an  instability  with  a  consequent 
large  increase  in  temperature  and  output  amplitude,  and  a  possible  propellant 
failure.  This  analysis  and  the  calculated  response  c  .rves  were  given  in  the 
previous  quarterly  report. 
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Experimental  data  confirming  the  existence  of  this  response  instability  were 
obtained  during  the  present  reporting  period.  Typical  experimental  results 
are  presented  in  Figures  4-7,  4-8,  and ,4-9.  The  data  represent  response 
parameters  at  each  frequency  immediately  after  vibration  is  started  and 
before  significant  heating  has  occurred  (isothermal  response)  and  also  the 
response  after  sufficient  time  to  attain  thermal  equilibrium  in  the  shear  speci¬ 
men  with  isothermal  boundaries  (equilibrium  reoponse).  The  theoretical 
responsi*  curves  for  the  steady- state  amplitude  ratio  and  insulated  surface 
temperature  are  reprinted  in  Figures  4-10  and  4-11.  from  the  first  quarterly 
technical  report  (Ref.  1)  for  qualitative  comparison.  The  theoretical  curves  - 
predict  response  for  a  hypothetical  propellant  with  an  isothermal  complex 
^  shear  modulus  log-log  slope  of  0.30.  Quantitative  comparison  of  the  experi¬ 
mental  data  with  theoretical  prediction  based  oxi  measured  propellant  proper¬ 
ties  has  been  completed  (with  reasonably  good  agreement  found)  and  vdll  be 
presented  in  a  subsequent  report. 


Interpolation  between  the  experimental  data  points  presented  in  Figures  4-7 
through  4-9  are  based  upon  the  theoretical  predictions.  The  experimental 
frequencies  which  can  be  used  are  limited  by  the  belt-gear  drive  mechanism 
of  the  milling  machine  used  for  the  vibration  tests. 

Of  particular  interest  is  the  double -valued  response  observed  in  the  steady- 
state  equilibrium  response  at  660  rpm  for  the  0.050 -inch  peak  input  displace¬ 
ment.  It  should  be  emphasized  that  each  equilibrium  data  point  represents 
the  sample  response  or  temperature  after  equilibrium  has  been  attained  at  the 
particular  vibration  frequency.  This  vibration  time  can  vary  from  20  minutes 
to  over  one  hour  for  the  teat  configuration  used.  For  example,  the  test 
sequence  for  the  0.050-inch  peak  displacement,  equilibrium  tests  was  as 
follows. 


With  the  sample  initially  at  room  ambient  temperature,  vibration  was  begun 
at  210  rpm  and  continued  until  thermal  (and  mechanical  response)  equilibrium 
was  attained.  The  excitation  frequency  was  then  immediately  raised  to  325 
rpm  and  maintained  until  a  new  equilibrium  was  established  and  so  forth, 
monotonically  increasing  the  excitation  frequency  in  steps  and  allowing  the 
specimen  time  tj  reach  equilibrium  at  each  frequency.  Alter  the  maximurr. 
frequency  equilibrium,  the  frequency  was  monotonically  decreased  in  steps, 
again  allowing  equilibrium  to  be  attained  at  each  frequency.  During  the 
increasing  frequency  sequence,  the  response  parameters  reached  equilibrium 
at  the  lower  660  rpm  data  point  for  amplitude  ratio,  insulated  surface  tem¬ 
perature  and  phase  angle.  However,  the  system  attained  equilibrium  at  the 
upper  branch  point  for  the  decreasing  frequency  sequence.  Prolonged  equi¬ 
librium  viuration  at  this  upper  branch  point  resulted  in  eventual  specimen 
failure  in  the  form  of  a  crack  propagating  along  the  direction  of  shear  and 
at  the  insulated  surface.  Inspection  of  the  failure  indicated  that  the  crack 
initiated  in  the  specimen  center  (the  highest  temperature  region)  and  propa  ¬ 
gated  outward  toward  the  semicircular  end  fillets. 
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V,  =  input  sinusoidal  DiSPLACtrf/LNT  AWPLITUDL 
Vm  =  OUTPUT  SINUSOIDAL  DISPLACEMENT  AMPLITUDE 
OF  MASS 

w  =  FREQUENCY 

wn  =  NATURAL  FREQUENCY  AT  AMOIENT  1  Cf.'.rERATURE 
V,/H=  NON-DIMENSIONAL  INPUT  DISPLACEMENT 
n  =  LOG-LOG  SLOPE  OF  ISOTHERMAL  COMPLEX  SMEAR 


MODULUS 


a>n 


Figure  4-10  Steady-State  Displacement  Amplitude  Ratio 
for  n  =  C.30 
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As  a  practical  result  of  these  investi^atious,  it  is  cle.ir  tliat  rocket  motor 
vibration  response  and.  anal/.sis  will  be  severely  coniplicated  by  thermo¬ 
mechanical  coupling  effects  if  the  inotor  configuration  and  \dbration  environ¬ 
ment  result  in  appreciable  inertial  loading  of  the  propellant  grain  as  a  whole, 
or  in  localized  areas  within  the  grain  (such  as  at  the  bar4e  of  a  slender  star 
point).  In  addition  to  the  usual  frequency,  a,mbient  temperature  and  vibration 
time  parameters,  it  is  clear  that  the  thermal  and  vibration  environncntal 
history  of  the  motor  must  be  considered  in  order  to  specify  and  define  the 
motor  requirements  and  vibration  response  characteristics.  It  would  seem 
that  the  length  of  time  which  must  be  considered  in  predicting  response  will 
be  of  the  same  order  as  the  thermal  time  constant  in  the  critical  region 
within  the  motor  grain. 

The  experimental  and  theoretical  results  show  the  multiple -valued  equilibrium 
response  region  to  occur  in  a  very  narrow  frequency  range  which  is  dependent 
on  the  vibration  loading  level,  the  natural  frequency  of  the  system  at  the  refer¬ 
ence  temperatu»'e ,  the  thermal  boundary  conditions,  and  the  propellant  visco¬ 
elastic  properties.  For  a  given  vibration  loading  level,  the  most  stringent 
loading  conditions  are  applied  to  the  propellant  along  the  upper  branch  of  the 
instability  response  region.  This  is  where  the  amplitude  ratio  is  a  maximum 
and  the  phase  angle  is  approximately  90  degrees,  resulting  in  large  strain 
amplitudes  in  the  propellant  as  well  as  high  equilibrium  temperatures. 

At  frequencies  above  this  narrow  region,  the  amplitude  ratio  decreases 
rapidly,  resulting  in  lower  strains  even  though  the  phase  difference  continues 
to  increase.  Below  the  instability  region,  propellant  strain  and  temperature 
are  both  low.  Both  lap-shear  propellant  specimens  used  for  the  investigation 
of  the  jump-instability  respon.se  failed  on  the  uppe"  branch  of  the  0.05-inch 
peak  input  displacement  instability  region  a*  660  rpin,  as  described  previously. 
For  the  one-inch  shear  specimen  gage  length  and  roughlv  90-degree  phase 
shift  at  that  point,  the  specimen  was  subjected  to  peak  strains  between  5  and 
10  percent,  consistent  with  vibration  failure  strain  levels  observed  previously 
for  constant  strain  amplitude  loading  conditions  (Ref.  1  ). 


4.3.3  Transient  Ther moviscoelastic  Response  for  Constant  Cyclic  Strain 
Amplitude  and  Inertial  Loading 


The  preceding  discussion  of  equilibrium  cyclic  inertial  loading  considered 
only  the  material  response  immediately  after  the  initiation  of  vibration 
(isothermal  response)  and  at  thermal  mechanical  steady-state  response 
(equilibrium  response).  Transient  thermal  and  mechanical  response  of  the 
inertially -loaded  system  is  very  difficult  to  calculate  analytically,  particu¬ 
larly  in  the  region  of  instability  where  the  mechanical  response  sweeps 
through  the  resonance  region  as  the  temperature  changes. 

To  investigate  the  effect  of  various  leading  and  thermal  boundary  conditions 
on  this  transient  rr  sponse,  a  series  of  experimental  tests  were  performed 
on  samples  identical  to  those  used  for  the  previous  equilibrium  tests. 
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Two  double -lap  shear  propellant  specimens  were  prepared  from  PBAA  pro¬ 
pellant  bars  milled  from  the  same  bulk  casting  (P2916).  One  sample  was 
prepared  with  wooden  support  members  and  compl€;tely  insulated  to  approxi¬ 
mate  adiabatic  conditions;  the  other  was  prepared  with  liquid-cooled  metal 
side  supports  and  a  low -conductivity  cent'^r  support  to  approximate  uni¬ 
directional  heat  flow,  Blach  specimen  was  subjected  to  a  series  of  four 
vibration  i  uns  --  tv/o  runs  with  constant  strain  amplitude  loading  at  2720  rprn 
and  2.47  percent  peak  sti-ain,  and  then  two  runs  with  inertial  loading  at  1115 
rpm  and  0.0247-inch  peak  input  displacement  amplitude.  The  multiple  ruriS 
were  performed  to  monitor  degradation  effects  due  to  the  combined  thermal- 
mechanical  environment.  Material  prooerties  were  measured  before  and 
after  each  vibration  run  at  a  frequency  3f  80  rpm,  room  ambient  tempter,  turc 
and  2.48  percent  peak  strain  amplitude.  The  results  of  these  single  frequency 
check  rums,  plotted  as  a  function  of  total  vibration  cycles,  are  presented  in 
Figvires  4-12  and  4-13  for  the  adiabatic  and  isothermal  boundary  conditions, 
respectively.  Appreciable  propellant  degradation  was  obtained  during  the 
vibration  sequence  as  evidenced  by  decreasing  storage  and  los.s  moduli.  The 
rate  of  degradation  is  seen  to  increase  during  the  inertial  loading  sequence 
due  to  the  larger  transient  strain  amplitudes  experienced  by  the  samples. 

Figures  4-14  and  4- 15  present  the  transient  mechanical  and  'hermal  response 
for  vibrations  runs  2  and  4,  referred  to  in  Figure  4-12,  for  the  adiabatic 
boundary  conditions.  Similarly,  Figures  4-16  and  4-17  present  the  same 
information  for  vibration  runs  2  and  4  referred  to  in  F'igure  4-13  for  the 
isothermal  boundary  conditions. 

i'he  storage  modulus,  loss  modulus  a.nd  loss  tangent  data  presen'  id  here  were 
calculated  u.sing  an  average  phase  angle,  due  to  the  nonlinearity  of  the  Lissajous 
pattern  ellipse  discussed  in  the  first  quarterly  report  (Ref,  1  ,  page  4-24, 
Figure  4-16),  The  tangent  of  the  angle  as  calculated  from  the  vertical  (force) 
intercepts  was  averaged  with  the  tangent  calculated  from  the  horizontal 
(displacement)  intercepts.  This  average  value  apjiears  to  correspond  most 
closely  with  the  effective  loss  tangent  as  calculated  from  the  measured  area 
of  the  Lissajon<?  or  hy.stercsis  pattern. 
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Section  5 

SUBCONTRACT  EFFORT 


5.1  PURDUE  UNIVERSITY  SUBCONTRACT  (Qu.-irterly  Summary  Report  to 
Lockheed  Propulsion  Company,  September  El  -  December  15,  1964. 
Professor  li.  A  Schapery  -  In  Charge) 


5.1.1  Nonlinear  Viscoelastic  Behavior  of  Highly  Filled  Propellants 


The  first  phase  of  this  program  is  concerned  v;ith  the  development  of  realistic 
stress- strain  equations  which  include  the  effects  of  strain  Induced  anistropy, 
but  are  sufficiently  simi>lc  to  permit  ea.sy  application  to  idealized  motor 
problems  in  the  second  phase  of  the  research.  Work  during  the  first  quarter 
dealt  primarily  with  the  evaluation  of  existing  theories  of  nonlinear  constitu¬ 
tive  equations  for  use  in  solid  propellant  applications  Nonlinear  equations 
range  from  extremely  general  forms  for  viscoelastic  behavior*  to  piecewise 
linear  approximations  for  clastic  response  *  In  the  former  case,  a  consider¬ 
able  number  of  mathematically  (rather  than  physically)  motivated  assumptions 
must  be  made  to  put  the  equations  in  a  practically  useful  form.  On  the  other 
hand,  the  latter  equations  of  piecewise  linear  elasticity,  with  coefficients 
depending  on  the  state  of  stress  or  strain,  seem  to  be  useful  engineering 
approximations  for  elastic  behavior,  however,  it  ia  not  at  all  clear  how'  these 
equations  can  be  generalized  to  vi.scoclasticity ,  except  by  trial  and  error, 
along  with  considerable  experimentation. 

In  another  paper,*  the  role  ol  thermodynamics  in  nonlinear  viscoelasticity 
was  Investigated,  and  it  was  shown  that  thermodynamic  principles  can  be  used 
to  develop  relatively  general  (but  simple)  stress- strain  equations.  Although 
Ref.  3  gives  only  a  onc-di  ncnsional  illustration,  wc  have  recently  made  the 
generalization  to  three  dimensional  behavior,  the  thermodynamic  basis  per¬ 
mitted  a  straightforward  generalization.  Briefly,  this  theory  makes  use  of 
the  first  and  sec<jnd  ’aws  of  thermodynamics  and  Onsager'e  principle.  Also, 
it  is  assumed  that  the  material  behaves  thermodynamically  as  a  linear  visco¬ 
elastic  body,  but  with  nonlinear  measures  of  deformation  in  place  of  the 
class  ical  linear  strains  The  stress  -  strain  equation.^  so  developed  contain 
a  natural  gene-alization  of  the  Boltzman  superposition  principle.  Moreover, 
the  equations  appear  to  be  consistent  with  existing  propellant  data,  and,  for 
example,  yield  a  relaxation  rhodulus  which  has  the  commonly  observed  form 
in  which  strain  and  time  dependence  appear  as  separate  factors. 

. . .  . . 

See,  for  example,  A.  IJ  Green  and  R  S.  Rivlin,  Archive  for  Rational 
Mechanics  and  Analysis,  \ol.  1,  pp  1-2*.  1957. 

^S.  B.  Dong,  L.  R  Herrmann,  K  S  Pister.  R.  L  Taylor,  Studies  Relating 
to  Structural  Analy.sis  of  Solid  Propellants,  Institute  of  Engineering  Research, 
"University  of  Califcirnia,  Berkeley.  1962'  ” 

^R.  A.  Schapery,  Journal  of  Applied  Physics,  Vol.  35,  No.  5,  pp  1451-  1465, 
1964. 
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Current  efforts  are  directed  at  applicatic)n  of  the  equations  to  the  mechanical 
property  test  geometries  used  by  1,PC;  in  order  to  evaluate  the  relatively  few 
unknf)v.'n  functions  and  to  determine  tl.e  limits  of  validity  of  the  theory. 

S.1.2  Thermomecl'.anical  Response  Studies  of  Solid  Propellants  Subjected  to 
Cyclic  and  Random  Loading 

A  report  )n  the  f i r s t -  qua r ter  work  has  heen  prepared  jointly  with  Dalton 
Cantey  of  LPC,  and  t\r  jicesentat  i!.>n  at  the  AlAA  6th  Solid  Propellant  Rocket 
Conference  to  be  held  February  1-3,  196S.  Briefly,  it  contains  an  extension 
of  [irevious  theoretical  studies  on  heat  gt'neration  to  ranrlom  loading,  provides 
solutions  to  the  thermal  and  mechanical  response  of  specimens  being  tested 
by  LPC,  compares  theory  with  experimental  data  (witli  reasonably  good  agree¬ 
ment  found),  and  draw.s  certain  praclii  al  conclusions  concerning  the  engineer¬ 
ing  usefulness  of  the  two  limit  cases  of  adiabatic  and  steady- state  thermal 
beha  V  io  r . 

Future  theoretical  investigations  will  study  the  the rmoic-echanical  behavior 
of  a  cantilever  beam  witli  and  v.'ithout  inertial  loading  Also,  we  shall  attempt 
to  improve  the  method  used  tc.»  preebet  heat  generation  by  taking  into  account 
nonlinearities  in  the  stress  -  strain  law. 
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5.2  UNIVERSITY  OF  WASHINGTON 


The  following  outlines  the  objectives  of  the  work  to  be  performed  at  the 
University  of  Washington  under  the  directi  n  of  Professors  R.  J.II.  Ballard 
and  E.H.  DHl. 

(1)  To  obtain  quantitative  data,  from  direet  observation  in  the  photo - 
viscoela-tic  bench,  of  the  effect  of  local  stress  concentrations  under  dynamic 
loadings  with  small  strains  and  constant  strain  rate  tests. 

(2)  To  investigate  experimentally  some  of  the  facets  of  the  nonlinear 
theory  for  viscoelastic  materials  which  has  been  proposed  by  Coleman  and 
Noll  (Foundations  of  Linear  Viscoelasticity,  Rev.  Mod.  Phys.,\^ol.  33(1961), 
pp.  239-249)  and  investigated  further  by  G.Lianis  in  previous  LPC  work. 

With  consideration  of  incompressible  fluids,  it  has  been  hypothesized  that 
some  of  the  material  constants  employed  in  the  nonlinear  constitutive  relations 
(large  strains)  can  be  determined  either  from  simple  dynamic  modulii  tests 

at  small  strains  or  creeptests  to  large  strains.  If  the  same  result  can  be 
shown  to  be  true  in  the  case  of  viscoelastic  solids  it  would  be  of  great  im¬ 
portance  to  the  definition  of  future  testing  programs  for  the  characterization 
of  solid  propellant  materials  in  the  nonlinear  range  of  deformation.  Since 
the  same  line  of  reasoning  which  led  to  the  hypothesis  concerning  the  constitu¬ 
tive  relation  can  be  applied  to  the  relation  of  the  dielectric  tensor  to  strain, 
this  relation  could  be  checked  by  simple  photoviscocTastic  observation.  This 
experimental  verification  of  the  above  stated  hypothesis  will  be  undertaken, 
employing  the  optical  bench  and  materials  now  available. 

Progress  during  the  reporting  period  was  as  follows.  The  photoelastic  bench 
has  been  modified  to  accept  a  double-lap  shear  specimen  of  viscoelastic  bire- 
fringent  plastic  and  to  provide  variable  frequency  sinusoidal  strain  inputs  of 
varying  amplitude-  Qualitative  observations  of  the  effect  of  the  stress  con¬ 
centrations  are  being  made. 

Theoretical  studies  of  viscoelastic  solids  have  been  carried  out  and  have 
raised  the  question  as  to  whether  or  not  the  nonlinear  constitutive  laws  can 
be  obtained  from  either  simple  dynamic  moduli  tests  at  small  strains  or 
from  creep  tests  to  large  strains.  These  theoretical  studies  are  near  com¬ 
pletion  and  will  be  reported  in  the  near  future. 
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Section  o 

FUTURE  PROGRAM 

Terhni  al  objectives  of  the  third  quarter  program  arc  listed  below. 

•  Completion  of  the  viscoelastic  response  and  multiaxial  failure 
investigation  of  the  remaining  group  of  high  solids  forn)Lilation 
variations  (80  percent  complete). 

•  Extension  of  multiaxial  failure  testing  and  investigation  of 
propellant  failure  surfaces  in  stress  space  (30  percent  complete). 

•  Deternnination  of  bulk  moduli  of  propellant  conspos itions  as  a 
function  of  applied  hydrostatic  pressure  and  temperature  (in 
progress) . 

•  Experimental  investigation  of  nonlinearity  effects  in  propellants 
under  various  types  of  test  for  development  of  realistic  non¬ 
linear  St  res s - strain  equations 
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